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Floquet-Engineered Valleytronics in Dirac Systems 
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Valley degrees of freedom offer a potential resource for quantum information processing if they can 
be effectively controlled. We discuss an optical approach to this problem in which intense light breaks 
electronic symmetries of a two-dimensional Dirac material. The resulting quasienergy structures 
may then differ for different valleys, so that the Floquet physics of the system can be exploited 
to produce highly polarized valley currents. This physics can be utilized to realize a valley valve 
whose behavior is determined optically. We propose a concrete way to achieve such valleytronics in 
graphene as well as in a simple model of an inversion-symmetry broken Dirac material. We study the 
effect numerically and demonstrate its robustness against moderate disorder and small deviations 
in optical parameters. 


Introduction .—Since the advent of graphene as a two- 
dimensional electronic material which can be produced 
in the laboratory BE], the possibility of exploiting the 
valley degree of freedom within it [3] and other Dirac 
systems has been vigorously studied. An important com¬ 
ponent of such valleytronic systems is the transport and 
detection of valley currents. Many of the ideas proposed 
to do so involve controlling the structure of the system, 
either using specific edge structures [3] or bulk nanos¬ 
tructures 0|5]. These ideas are limited by the precision 
they require to control the nanostructure. A particularly 
interesting way around these limitations combines intrin¬ 
sic band properties with optics to yield valley-contrasting 
behavior. In gapped Dirac systems such as M 0 S 2 and 
WS 2 , or bilayer graphene in a perpendicular electric field, 
valley currents can be induced using the differing Berry’s 
curvatures of the valleys 1303 ]. In such systems, cir¬ 
cularly polarized light can excite different electron-hole 
pair populations in different valleys mm, leading to a 
“valley Hall effect” 0 9] [15HT71 . 

In this work we discuss a fundamentally different ap¬ 
proach to optically-controlled valleytronics in the bulk 
that offers a high degree of tunability in a single sample. 
In this approach, the light is relatively intense, so that 
the electronic structure is represented by eigenvalues of a 
Floquet Hamiltonian. The time-dependence of the elec¬ 
tric field, rather than intrinsic properties of the material 
or nanoscale structures, is used to effectively break in¬ 
version symmetry and distinguish the valleys. One way 
to do this, for example in graphene, is by shining an 
admixture of circularly polarized light of frequencies Ul 
and 4$2, which can be coherently generated by use of 
nonlinear crystals. As explained below, for appropriate 
choices of amplitudes and phase offset one can produce 
a Floquet quasienergy spectrum which is gapped for one 
valley but gapless for the other. A dc current passed be¬ 
tween leads with chemical potentials in this gap is then 
valley-polarized. The degree of polarization can be in¬ 
terrogated with different phase-offsets in the vicinities of 
each lead, such that the gap closing is in the same or dif¬ 
ferent valleys for each. An example of this behavior for 
an idealized system is illustrated in Fig. [lj The polar¬ 
ization turns out to be quite robust against disorder and 



FIG. 1. (color online) (a) The schematics of the graphene 
system: A background laser with frequency Q and intensity 
(c/ 47 r)AfH 2 is supplemented with a circularly-polarized 4th 
harmonic of frequency 4Q and intensity ( 4 c/ 7 t)A|D 2 with a 
relative phase offset 4 >l on the left and (j)R right half of the 
device, (b) The valley polarization Pq — 1 — cr 0 ff/cr on in 
a ribbon with periodic boundary conditions vs. length L 
(in units of the lattice constant no)- The inset shows the 
on- and off-conductances, a Q ff and a on , respectively. Here, 
eAiao/c = 1.0, A 4 /A 1 — 0.1, n /7 « 1.47, where 7 « 2.7 eV is 
the hopping parameter, and the ribbon width w = 6^/3ao. 


edge effects, as shown below. This system represents an 
optically controlled valley valve. 

This physics can also be applied to systems in which in¬ 
version symmetry is already broken, such as graphene on 
a BN substrate [9] or dichalcogenide materials like M 0 S 2 
and WS 2 , which have pre-existing gaps that are the same 
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FIG. 2. (color online) (a) The bulk quasienergy spectrum in 
the Brillouin-Floquet zone for uniformly irradiated graphene 
(eAiao/c = 1.5, A 4 /A 1 = 0.6,0 = 7 t/ 2 ,Q /7 = 1.44) shows 
a gap at K_ and gapless Dirac dispersion at K + point, (b) 
and (c) The quasienergy gaps ClS± of the K + (blue) and K_ 
(orange) valleys at £ = Cl/2 show the gaps may be opened 
and closed by tuning the optical parameters. In (b) a\ = 
eAiao/c = 1.5, Cl/'y = 1 and in (c) 0 = 0 , Cl/'y = 1.44. 


for both valleys. Circularly polarized light may close one 
gap while opening the other in such materials, again al¬ 
lowing marked preferential conduction for one of the two 
valleys. Left- and right-circularly polarized light create 
open channels for opposite valleys, leading to optically- 
controlled valley polarization, as we demonstrate below. 

Optically Broken Inversion Symmetry in Graphene — 
In the presence of a temporally periodic potential, 
electronic states follow the time-dependent Schrodinger 
equation. Floquet’s theorem m guarantees that its so¬ 
lutions as a function of time t have the form ifj a (t) = 
u a (t)e~ l£cxt , with u a (t + T) = u a (t ), where T is the pe¬ 
riod of the Hamiltonian H(t), and a includes any quan¬ 
tum numbers required to specify the electronic state. The 
quasienergies (which may be restricted to the inter¬ 
val — Cl/2 < e a < Cl/2) are eigenvalues of the “Floquet 
Hamiltonian” Hp(t) = H{t) — id t , and u a are the corre¬ 
sponding eigenfunctions. 

As a paradigm of these systems we consider the case 
of graphene, in which the electrons are assumed to hop 
around on a tight-binding honeycomb lattice [2lfl9l [20] . 
and the effects of the circularly polarized light are imple¬ 
mented via time-dependent phases in the hopping ma¬ 


trix elements [ 21 ]. Variations in the intensity and the 
frequency of the light tune the quasienergy band struc¬ 
ture through many distinct topological phases [22] . The 
energy spectrum includes two valleys of states near the 
K + and K_ points of the Brillouin zone, and wavefunc- 
tions are two-component spinors representing the elec¬ 
tron amplitudes on the two sublattices of the honeycomb 
structure. In the presence of a spatially uniform, time- 
dependent electric field from light normally incident on 
the graphene plane, the Floquet Hamiltonian in the sub¬ 
lattice basis {ua,ub) has the form [ 22 ] 


H F (k,t) 


-id t — 7 Z(k,t) \ 

-7 -id t )' 


( 1 ) 


where k is the wavevector of the state, 7 is a hopping am¬ 
plitude, Z(k, t) = Y^n= 1 eh k +f A ( t )]' ari 5 a n are the near¬ 
est neighbor vectors of a site on the lattice, and A(t) = 
Ai (cos Clt, sin Clt) is the vector potential for the electric 
field, with Cl = 2i r/T. When k is set to the K± point 
[K± = (0, ±4:7r/3\/3ao ), with ao the nearest neighbor dis¬ 
tance] the lattice symmetry combines with the temporal 
symmetry so that Z(K±,t + T/3) = e^ 2 ™/ 3 Z(K±, t), ef¬ 
fectively tripling the relevant frequency for the Floquet 
problem. Because of this, distinct states can cross, rather 
than repel, at the Floquet zone boundary = ±Cl/ 2 as 
A\ or Cl are varied, leading to topological transitions in 
the quasienergy band structure. 

In this situation such crossings occur at the K± points 
simultaneously. This results from a combination of inver¬ 
sion symmetry (K_ = — K+) and a form of time reversal 
symmetry: If u(K+,t) is an eigenvector of Hp(K+,t), 
then a x u*( K + , —t), with cr x a Pauli matrix, is an eigen¬ 
vector of Hp(K-,t) with the same eigenvalue. Lifting 
this coincidence of eigenvalues distinguishes the valleys. 
One way to do this is by changing Hp such that 


Z(K+,-i)^Z(K_,t), (2) 

while retaining the T/3 period of e ±int Z(K±,t) needed 
for a gap closing, breaking the effective inversion symme¬ 
try. In contrast to approaches which use static potentials 
SIS], here we seek to use the vector potential to do so, 
allowing for optical control of the valley-distinguishing 
properties of the system. 

One way this can be done is by adding a 4th har¬ 
monic to the vector potential, so that it has the form 
A x +iA y = Aie lClt + A. 4 e^ 4m+< ^. For this A the inequal¬ 
ity ( 2 ]) is satisfied, provided 0 ^ 0 , 7 r, while the frequency 
tripling of e ±1,nt Z(K±,t) is retained. As illustrated in 
Fig. 11 gaps at K± are in general unequal and can be 
opened and closed around e = Cl/2 separately by tuning 
the optical parameters. Note that the phase offset 0 may 
be adjusted by varying the optical path length of the 4 Cl 
light component relative to the Cl component. As we next 
show numerically, this leads to different band gaps for the 
two valleys, even allowing the gap to close for one while 
the other remains open. Thus, the admixture of the two 
frequencies of light in principle allows one to prohibit a 
bulk current for one valley while allowing it for the other. 
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Numerical Results .—To test this idea, we have com¬ 
puted the DC conductance at zero temperature for an 
irradiated graphene strip [22, 23;. Armchair graphene 
ribbons were simulated, with periodic boundary condi¬ 
tions across the width to minimize the edge effects while 
keeping the size of the system small enough for numer¬ 
ical efficiency. Optical parameters were chosen to pro¬ 
duce quasienergy gaps large enough to observe the effect 
with our simulated system sizes. More realistic optical 
parameters are discussed below. We also performed sim¬ 
ulations for more realistic systems with open boundary 
conditions. These results are reported below and support 
our conclusions (see Appenix ©• 

The current is introduced from the leads attached to 
the two ends of the system, while the central graphene 
region is voltage-biased to align the Floquet zone edge 
with the average chemical potential of the leads. In this 
scheme, the entire system is illuminated by light of fre¬ 
quency Q. Each half of the system around the leads is fur¬ 
ther illuminated by circularly-polarized light of frequency 
4f] with independently controlled values of the phase off¬ 
set, (J)l and (j)R . With </>£ = </>#, due to the gap for one 
valley, the system essentially allows only current from the 
other valley to pass, yielding a valley-polarized current. 
That this current is valley-polarized is confirmed by com¬ 
paring with the conductance when (f>L — &r = (2 n + 1)7r, 
with n an integer. In this case, the two halves are conduc¬ 
tive for opposite valleys so very little net current passes. 
This means the system behaves as a valley valve [3. which 
may be opened or closed optically. The relative conduc¬ 
tance in these two cases (cr on and cr 0 ff, respectively) offer 
a measure of the valley polarization, P 0 = 1 — <j 0 &/(j 0n , 
achieved in the system. The fidelity of the resulting valve 
can exceed Pq = 98%, as illustrated in Fig. [ljb) (see also 
Appendix [D] This is one of our main results. 

Effects of Edges .—The valley polarization in our pro¬ 
posal results from bulk transport. Due to the existence 
of multiple Floquet topological phase transitions tuned 
by frequency [22], the system with open boundary condi¬ 
tions has, in addition, a number of chiral edge states. 
These edge states can carry current but lack a well- 
defined valley index. Moreover, they can scatter the bulk 
states between valleys. Both of these effects can degrade 
the observed valley polarization, and are particularly no¬ 
ticeable in small systems. We expect that in sufficiently 
large sample the degrading effects of the edges will be 
much reduced. 

To support this expectation, we also simulated a more 
realistic system with open boundary conditions. To min¬ 
imize edge effects, we employ a geometry in which the 
leads are connected away from them, as illustrated in 
Fig. 0 a). This requires relatively large widths, limiting 
the system lengths one can ultimately study efficiently. 
Nevertheless, one may still obtain information about the 
large-length limit by choosing a quasienergy gap that 
is not too small. We report our results in Fig. [3] and 
compare to the case with periodic boundary conditions. 
The path followed in Fig. [3|a) is chosen so that the gap 
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FIG. 3. (color online) The valley polarization Po for the sys¬ 
tem with periodic (dotted blue line) and open (solid green 
line) boundary conditions, respectively. The index i indicates 
parameter values at points along the $+ = 0 (orange) contour 
shown in the inset, as in Fig. [2^c) . The corresponding gap 
5- is also shown (dashed black line). The fixed parameters 
are as in Fig. [2ja) ; the chemical potential of the irradiated 
region is 0.717 ~ Q/2; the periodic system has a full length 
L = 2£ = 288ao and width w = 6V3ao] the open system has 
£ = 27ao, w = 24v / 3«o; the leads have a width 12v / 3«o and 
are connected equidistant from the edges across the width and 
a distance 9ao away from the center. 



FIG. 4. The polarization of the graphene valley valve vs. the 
strength of onsite static disorder potential, W. The param¬ 
eters are as in Fig. [2ja) and the chemical potential of the 
irradiated region is 0.7Fy ~ Q/2. The boundary conditions 
are open and the dimensions are as for the open system in 

Fig-i 

for one of the valleys vanishes precisely, maximizing the 
“on” current and thereby the valley polarization. The 
system with periodic boundary conditions and a gapless 
K + point shows nearly 100% valley polarization for an 
arbitrary value of the quasienergy gap L > 0 at K_. 
For the system with open boundary conditions, a larger 
gap 5- is required to have significant valley polarization. 
Even so, remarkably large valley polarizations are ob¬ 
tained in this case for a relatively small system (see also 
Appendix [B| . 

Dirac Systems with Statically Broken Inversion Sym- 
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metry. —A second way to break the symmetry of the val¬ 
ley electronic states and realize an optically-controlled 
valleytronic system is by statically lifting the inversion 
symmetry of the system. Examples of such systems are 
provided by graphene deposited on boron nitride [24] [25] 
and single-layer M 0 S 2 mug. a simple description of 
these systems is a tight-binding model on a honeycomb 
lattice with a staggered potential d= fi s of opposite signs 
on different sublattices; the low energy electron states 
at the K± points are then governed by massive Dirac 
Hamiltonians. Despite this apparent inversion symme¬ 
try breaking, the magnitude of the gap in both valleys 
is the same and no valley polarization can be realized in 
equilibrium. Valley polarization in a static region can be 
achieved in this system via optical absorption P31EZI- 
In our approach, the gaps at the K± points can 
be distinguished by the helicity of the monochromatic 
circularly-polarized light. The Floquet Hamiltonian for 
the low-energy excitations now takes the form Hp( k, t) + 
M, where Hp is the same as in Eq. ([!]), and 

M =(o- 0 ,.)- < 3) 

Due to the inversion-symmetry breaking, the quasiener¬ 
gies around the two valleys now evolve differently: for a 
given helicity, the frequency can be tuned to a value such 
that there is a gapless K + valley crossing 6 = 0, and a 
gapped K_ valley (see Appendix 0- Illuminating the 
two halves of the system with circularly-polarized lasers 
at such a frequency realizes a valley valve. The valve can 
be turned on and off by switching the helicity of one of 
the lasers. We have checked that, in this case, the cur¬ 
rent induced from the leads at energy E « 0 has valley 
polarization exceeding 90%. 

Discussion. —To observe the effects we have demon¬ 
strated above, energy and length scales in real samples 
must be chosen appropriately. For example, the temper¬ 
ature in the leads must be less than the quasienergy gap 
of the gapped valley, namely, Dc)_. This can be tuned op¬ 
tically. For example, at the main frequency D = 0.2y = 
0.54 eV, a gapless K + is obtained at a% = eA^a^fc = 0.1 
and laser intensity I\ = ( 1 / 47 xc)A\D? ~ 10 14 W/m 2 . At 
these values, by varying the ratio A 4 /A 1 , the gap at K_ 
can be tuned to be > 5 meV. Increasing the frequency 
and the intensity of the laser can produce even larger 
gaps (see Appendix 0 . Since the conductance of the 
gapped valley decays over the wavefunction evanescent 
length, £ ev oc v/Q5-, the length of the system L must 
exceed £ ev . Here, v is the Fermi velocity at the valley, 
which can also be tuned optically. In our estimate, for 
the aforementioned values, £ ev < 500ao = 71 nm 

In practice, there are always extrinsic perturbations 
limiting the polarization. Prominent among these is dis¬ 
order, which does so through inter-valley scattering. In¬ 
tuitively, the effect of disorder is to fill in the gap in 
the quasienergy spectrum. We expect that its effect is 
controlled by this gap, and that it should not spoil the 
valley polarization when weak enough. In Fig.[4]we show 


results of our simulations of disorder-averaged polariza¬ 
tions in the graphene system discussed above (see also 
Appendix JC|. In these simulations we have used ran¬ 
dom, static onsite potentials, with Gaussian distribution 
of zero mean and standard deviation W characterizing 
the disorder strength. The conductances are each aver¬ 
aged over 25 disorder configurations. Indeed, at larger 
strengths, disorder causes intervalley scattering and di¬ 
minishes the valley polarization by increasing cr 0 ff (see 
inset). However, when W is sufficiently smaller than the 
quasienergy gap, here Dc)_ ~ O.I 7 , the valley polariza¬ 
tion we obtain is robust against disorder, exceeding 70% 
over a significant range of disorder strength. As discussed 
above, in real samples the quasienergy gap can be tuned 
to be > 5 meV. Samples with disorder strengths below 
this appear to be currently available in the lab pTH34] . 
Thus we believe there is considerable room to vary the pa¬ 
rameters of the system without having the effects spoiled 
by disorder. 

Novel forms of electronic states arise as the result of the 
interaction between electrons and periodic external driv¬ 
ing fields. The optically controlled valley polarization 
and valleytronic devices proposed in this work promise 
new ways of engineering and utilizing the emergent elec¬ 
tronic degrees of freedom in graphene and related Dirac 
systems. 

This work was supported in part by the NSF through 
Grant Nos. DMR-1350663 and DMR-1506460, the US- 
Israel Binational Science Foundation, and by Indiana 
University. 


Appendix A: Optical parameters 

Here we first show numerical results for the optical pa¬ 
rameters needed for the graphene valley valve discussed 
in the main text and in Fig. 1. In Fig. |5|a), we plot the 
solutions to 5 + (Ai,Ai,fi) = 0 with the largest D, i.e. 
the points at which the quasienergy gap at 6 = D /2 first 
vanishes for the K + valley. In Fig.[5|b) we plot the value 
of the gap at these points for K_ point. Together, these 
plots can be used to guide the choice of optical param¬ 
eters in an experiment. For example, laser intensity of 
the main frequency I\ = ( 1 / 47 tc)A 2 D 2 ~ 10 14 W/m 2 for 
oq = eAido/c = 0.1 and D = 0 . 2 q = 0.54 eV, which is 
within current capabilities. By adjusting the intensity of 
the 4th harmonic, one can tune the gap for the K_ point 
shown in Fig. [5|b) by more than an order of magnitude; 
for an intensity I 4 ~ I\ in graphene (7 « 2.7 eV), this 
gap is ~ 5 meV. 


Appendix B: Numerical methods 

We have employed the Green’s function method [28] 
for computing the conductance of the system numeri¬ 
cally. The recursive Green’s function method [29] is em¬ 
ployed for the computation of the density of states. The 
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FIG. 5. (a) The optical parameters a± = eAiao/c , A 4 /A 1 
and the largest fl /7 at which the quasienergy gap <5+ = 0 
at quasienergy e — Q/2 for the graphene device, (b) The 
corresponding gap at K_. 


leads were modeled as highly doped (to l/6th of their 
bandwidths) semi-infinite graphene ribbons. 

The system is described by the Hamiltonian Hs(t ) = 
H(t) + He , where H(t) = H(t + T) is the Hamiltonian 
of the irradiated area, and the contact Hamiltonian is 

He = ^ F x t a x ^ci + h.c. = ^ a x ^F x c J r h.c., 

A =L,R X=L,R 



FIG. 6. The polarization P(A<j>) = 1 — cr(7r)/cr(A0) for 
the graphene system with periodic (dotted blue line) and 
open (solid green line) boundary conditions, respectively, vs. 
A0 = (j)R — <j>L. Also shown is the conductance cr(A0) at 
zero temperature (dashed black line) for the system with open 
boundary conditions. The fixed parameters are as in Fig. 2(a) 
of the main text; the chemical potential of the irradiated re¬ 
gion is O. 7 I 7 ~ Q/2; the periodic system has a full length 
L = 2 £ = 288ao and width w = 6\/3&o; the open system has 
£ = 27ao, w — 24v / 3«o; the leads have a width 12v / 3«o and 
are connected equidistant from the edges across the width and 
a distance 9ao away from the center. 


where the Floquet Green’s function is 

G«( w ) = i J dtjds G(t,t — s)e iujs e ikQt , fee Z 

in terms of the time-dependent Green’s function of the 
system G(t, t') = G(t + T, t' + T). 

The coupling between the irradiated area and the 
leads enters through the factors £ a (cj) = 2ttF x ^p x (uo)F x , 
where p x ( cj) is the density of states in lead A. Finally, 
for zero relative bias of the leads, the zero-temperature 
conductance is given by 


where L, R denote left and right leads. The a x t and 
c\ denote, respectively, the electronic creation opera¬ 
tor at site r for lead A and at site l for the irradiated 
area. The time dependent current through lead A is 
J x (t ) = ie [H{t ) + He , A/' A (t)], with N x the number op¬ 
erator. The time averaged current in the left lead is 
/ = J Q T (J L (t))dt/T, where (•) denotes an average over 
the lead states. It may be written in the form [35, 

/ = h / dw E [tAh/^h - , 

where f\(uj) = [1+ e^ _eVA ^ TA ] _1 is the Fermi-Dirac dis¬ 
tribution in lead A, with electric bias V\ and temperature 
r\. (We adopt units in which the Boltzmann constant 
k,B = 1.) The transmission probability is given by 


(7 = £e[ t A(°)+Al( o)[ • 

The conductance computation is simplified by the use 
of the “wide band limit”, p x (uj) = p x independent of 
uj [22] . In Fig. [6j we show the results of the simulation 
for conductance and valley polarization as a function of 
the phase difference A</> = 4>r — <Pl, as in the setup shown 
in Fig. 1(a) of the main text, for both periodic and open 
boundary conditions. In either case, the system shows 
almost purely valley-polarized transport and a large “on” 
signal for a wide range around A</> = 0. 

Appendix C: Dirac system with broken inversion 
symmetry 


rp(F) _ 

1 XX' — 


Tr 


G (fc)t (w)f A (w + kn)G (k) (co) 


A schematic of the system to create valley-polarized 
current using statically broken Dirac system is shown in 
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FIG. 7. (a) Schematic device with two circularly polarized light of opposite helicity incident on a Dirac system with broken 
inversion symmetry (such as M 0 S 2 ). (b) The quasienergy spectrum for a uniformly irradiated system modeled with a staggered 
chemical potential with fixed helicity shows a gapless K + point at e — 0 and a gapped K_ point; The gapped and gapless points 
are switched for the opposite helicity. The staggered chemical potential fi s = O. 47 . The other parameters a± = eAiao/c = 
0.5, Q /7 = 1.73, and the chemical potential of the irradiated region is 0.1 27 . 



FIG. 8 . The polarization of a gapped Dirac system vs. the 
strength of onsite static disorder potential, W. The pa¬ 
rameters used are as in Fig. [9] but the system has open 
boundary conditions. The dimensions of the system are 
w = = 54ao, the leads have a width 7v / 3«o and 

connected a distance 36ao away from the center. 

Fig. Ma). A typical quasienergy spectrum is shown in 
Fig. Mb), where the inversion symmetry is broken by 
staggered potential, as in Eq. (3) of the main text. Also, 
in comparison with Fig. 4 of the main text, we show in 
Fig. [8] the effect of disorder on the valley polarization in 
this case. 


Appendix D: Length dependence of the valley 
polarization 

As seen in Fig. 1(b), the valley polarization, Pq = 
1 — (Joff/^on, is negative for a range of lengths before ap¬ 
proaching 1 for larger lengths. In this section, we present 
a model for this length dependence. 

The reason for this particular behavior lies in the inter¬ 
play of diffusive and evanescent transport in a finite size 
system. As a simple model, we shall assume that the two 
valleys behave as two independent conduction channels 



FIG. 9. Modeling the valley polarization vs. length L of 
the system shown in Fig. 1(a). The width w = 12ao, the 
evanescent length £ ev — 20ao and the interpolating function 
f(x) — tanhx. 


in each half of the device shown in Fig. 1(a). The gap 
at one of the valleys, say K + , vanishes at large system 
sizes (diffusive conduction channel) while the other one, 
say K_, remains finite (evanescent conduction channel). 
In the “off” configuration, the valley currents are carried 
in parallel within each half of the system, and the two 
halves are connected in series. In the “on” configuration, 
the valley currents are carried in parallel along the whole 
length of the system. At small system sizes, the evanes¬ 
cent length becomes of order L, the distance between the 
leads. In this regime, both channels are diffusive. As 
the system size increases, only one of the parallel chan¬ 
nels in the “on” configuration remains diffusive; in the 
“off” configuration, both parallel channels behave simi¬ 
larly. This causes a faster drop in the “on” conductance 
in intermediate length scales. At yet larger length scales, 
the “on” configurations remains diffusive in one channel, 
while the “off” current becomes completely evanescent 
and decays exponentially. 

To make this more precise, we take the conductance a 















for each valley to be 


satisfying f(x <C 1) = x and f(x^> 1) = 1. The “on : 
and “off” conductances are given by 
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V 


cr(l,Z) = J exp 



(Dl) 


Here, w is the width of the system, l is the length of the 
channel, and £ is a length scale inversely proportional to 
the valley’s quasienergy gap. For the K + valley, trans¬ 
port is diffusive and we take = L. For the K_ valley, 
there is a crossover from diffusive to evanescent trans¬ 
port; so, we take = £ ev f(L/£ ev ), where £ ev is the 
evanescent length and f{pc) is an interpolating function 


cron = cr(L,tf) + cr (A£-), 


^off — ^ 


1 

al] 2 , e +) 


(t(L/2,£_)J 


(D2) 

(D3) 


Our results are plotted in Fig. [9] for a representative set 
of parameters, showing the same behavior for Po(L) in 
this model as in Fig. 1(b). 
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